(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 10.4' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     23292,        604]
NotebookOptionsPosition[     22650,        582]
NotebookOutlinePosition[     22991,        597]
CellTagsIndexPosition[     22948,        594]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[{
 RowBox[{
  RowBox[{"Vesc", "=", "2380.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"grav", "=", "1.622"}], ";"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Impact", " ", "target", " ", "material", " ", "parameters"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"K1", " ", "=", " ", "0.2"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mu", " ", "=", " ", "0.55"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Ybar", "=", 
   RowBox[{"18.0", "*", 
    RowBox[{"10", "^", "6"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"rhot", "=", "2500"}], ";"}], 
  RowBox[{"(*", "2680.0", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"rhoi", "=", "2500"}], ";"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Impact", " ", "conditions"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"ai", "=", "4000.0"}], 
   RowBox[{"(*", 
    RowBox[{"7000.0", ";", "1100.0", ";", "1600.0", ";", "850.0"}], "*)"}], 
   ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"ai", "=", 
     RowBox[{"4", " ", "km"}]}], ",", " ", 
    RowBox[{"vi", "=", 
     RowBox[{"18.3", " ", 
      RowBox[{"km", "/", "s"}]}]}], ",", " ", 
    RowBox[{"thetai", "=", 
     RowBox[{"90", " ", 
      RowBox[{"deg", ":", " ", "copernicus"}]}]}]}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mi", "=", 
   RowBox[{"4.0", "*", 
    RowBox[{"Pi", "/", "3.0"}], "*", "rhoi", "*", 
    RowBox[{
     RowBox[{"(", "ai", ")"}], "^", "3"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vi", "=", "18300.0"}], 
  RowBox[{"(*", 
   RowBox[{"12000.0", ";", "18000.0", ";", "17800.0", ";", "20000.0"}], 
   "*)"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"thetai", "=", 
   RowBox[{"Pi", "/", "2.0"}]}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Pi", "/", "2"}], ";", 
    RowBox[{"Pi", "/", "4.0"}]}], "*)"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vi", "=", 
    RowBox[{"vi", "*", 
     RowBox[{"Sin", "[", "thetai", "]"}]}]}], ";"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Impact", " ", "crater"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pi2", "=", 
   RowBox[{"grav", "*", 
    RowBox[{"ai", "/", 
     RowBox[{"(", 
      RowBox[{"vi", "^", "2"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pi3", "=", 
   RowBox[{"Ybar", "/", 
    RowBox[{"(", 
     RowBox[{"rhot", "*", 
      RowBox[{"vi", "^", "2"}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pi4", "=", 
   RowBox[{"rhot", "/", "rhoi"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vols", "=", 
   RowBox[{"K1", "*", 
    RowBox[{"mi", "/", "rhot"}], "*", 
    RowBox[{
     RowBox[{"(", " ", 
      RowBox[{
       RowBox[{"pi2", "*", 
        RowBox[{
         RowBox[{"(", "pi4", ")"}], "^", 
         RowBox[{"(", 
          RowBox[{
           RowBox[{"-", "1.0"}], "/", "3.0"}], ")"}]}]}], " ", "+", " ", 
       RowBox[{
        RowBox[{"(", "pi3", ")"}], "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"mu", "+", "2.0"}], ")"}], "/", "2.0"}], ")"}]}]}], " ", 
      ")"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "3.0"}], "*", 
       RowBox[{"mu", "/", 
        RowBox[{"(", 
         RowBox[{"mu", "+", "2.0"}], ")"}]}]}], ")"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"volg", "=", 
   RowBox[{"K1", "*", 
    RowBox[{"mi", "/", "rhot"}], "*", 
    RowBox[{
     RowBox[{"(", " ", 
      RowBox[{"pi2", "*", 
       RowBox[{
        RowBox[{"(", "pi4", ")"}], "^", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"-", "1.0"}], "/", "3.0"}], ")"}]}]}], " ", ")"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "3.0"}], "*", 
       RowBox[{"mu", "/", 
        RowBox[{"(", 
         RowBox[{"mu", "+", "2.0"}], ")"}]}]}], ")"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Rs", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"3.0", "*", 
      RowBox[{"vols", "/", "Pi"}]}], ")"}], "^", 
    RowBox[{"(", 
     RowBox[{"1.0", "/", "3.0"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Rg", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{"3.0", "*", 
      RowBox[{"volg", "/", "Pi"}]}], ")"}], "^", 
    RowBox[{"(", 
     RowBox[{"1.0", "/", "3.0"}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Df", "=", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"2.0", "*", 
        RowBox[{"Rg", "/", "1000.0"}]}], ")"}], "^", "1.18"}], ")"}], "*", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"(", "20.2", ")"}], "^", 
      RowBox[{"(", 
       RowBox[{"-", "0.18"}], ")"}]}], ")"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Yt", "=", 
    RowBox[{"rhot", "*", 
     RowBox[{"(", 
      RowBox[{"vi", "^", "2"}], ")"}], "*", 
     RowBox[{
      RowBox[{"(", " ", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"grav", "*", 
          RowBox[{"ai", "/", 
           RowBox[{"(", 
            RowBox[{"vi", "^", "2"}], ")"}]}]}], ")"}], " ", "*", " ", 
        RowBox[{
         RowBox[{"(", 
          RowBox[{"rhoi", "/", "rhot"}], ")"}], "^", 
         RowBox[{"(", 
          RowBox[{"1.0", "/", "3.0"}], ")"}]}]}], " ", ")"}], "^", 
      RowBox[{"(", 
       RowBox[{"2.0", "/", 
        RowBox[{"(", 
         RowBox[{"mu", "+", "2.0"}], ")"}]}], ")"}]}]}]}], ";"}], 
  "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Ejection", " ", "velocity", " ", "formula", " ", "parameters"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Ctg", "=", "0.85"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Cvpg", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"Sqrt", "[", "2.0", "]"}], "/", "Ctg"}], " ", "*", 
    RowBox[{"(", 
     RowBox[{"mu", " ", "/", " ", 
      RowBox[{"(", 
       RowBox[{"mu", " ", "+", " ", "1.0"}], ")"}]}], ")"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Cvps", "=", 
   RowBox[{"Cvpg", "*", 
    RowBox[{"Sqrt", "[", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"rhot", "*", "grav", "*", "Rg"}], ")"}], "/", 
      RowBox[{"(", 
       RowBox[{"Ybar", "+", "Yt"}], ")"}]}], "]"}], "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"Rg", "/", "Rs"}], ")"}], "^", 
     RowBox[{"(", 
      RowBox[{"1.0", "/", "mu"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"fric", "=", 
   RowBox[{
    RowBox[{"Ybar", "/", "rhot"}], "*", 
    RowBox[{
     RowBox[{"(", "Cvps", ")"}], "^", "2"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"ge", "[", "r_", "]"}], "=", 
   RowBox[{"r", "*", "grav", " ", "*", " ", 
    RowBox[{
     RowBox[{"(", "Cvpg", ")"}], "^", "2"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"ve", "[", "r_", "]"}], "=", 
   RowBox[{"Cvpg", "*", 
    RowBox[{"Sqrt", "[", 
     RowBox[{"grav", "*", "Rg"}], "]"}], "*", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"r", "/", "Rg"}], ")"}], "^", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"-", "1.0"}], "/", "mu"}], ")"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vef", "[", "r_", "]"}], "=", 
   RowBox[{"Sqrt", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"ve", "[", "r", "]"}], "^", "2"}], " ", "-", " ", 
     RowBox[{"ge", "[", "r", "]"}], " ", "-", " ", "fric"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"l", "[", "r_", "]"}], "=", 
   RowBox[{"r", " ", "+", " ", 
    RowBox[{"2.0", " ", "*", " ", 
     RowBox[{"vef", "[", "r", "]"}], "*", 
     RowBox[{"Sin", "[", 
      RowBox[{"Pi", "/", "4.0"}], "]"}], " ", "*", " ", 
     RowBox[{"vef", "[", "r", "]"}], "*", 
     RowBox[{
      RowBox[{"Cos", "[", 
       RowBox[{"Pi", "/", "4.0"}], "]"}], "/", "grav"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"r0", "=", "ai"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"lg", " ", "=", " ", "400000.0"}], ";", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Df", "*", "1000.0", "*", "4"}], ";"}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{"pos", "=", 
   RowBox[{"FindRoot", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"l", "[", "r", "]"}], "\[Equal]", "lg"}], ",", 
     RowBox[{"{", 
      RowBox[{"r", ",", "r0"}], "}"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"pos1", "=", 
   RowBox[{"Re", "[", 
    RowBox[{"r", "/.", "pos"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"pos1_scale", "=", 
    RowBox[{"pos1", "/", "Rg"}]}], ";"}], "\[IndentingNewLine]", 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Ballisitc", " ", "flight", " ", "time"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"thetae", " ", "=", 
    RowBox[{"Pi", "/", "4.0"}]}], ";"}], " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"Timeballistic", " ", "[", "pos1", "]"}], "=", " ", 
    RowBox[{"2.0", " ", "*", " ", 
     RowBox[{"vef", "[", "pos1", "]"}], "*", 
     RowBox[{
      RowBox[{"Sin", "[", "thetae", "]"}], "/", "grav"}]}]}], ";"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Cooling", " ", "time", " ", 
    RowBox[{"estimatio", ":", " ", 
     RowBox[{"SI", " ", "unitn"}]}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"temperature", " ", "of", " ", 
    RowBox[{"environment", ":", " ", 
     RowBox[{"room", " ", "temperature", " ", "in", " ", "Kelvin"}]}]}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"TempEnv", "=", " ", "200.0"}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Temperature", " ", "of", " ", "molten", " ", 
    RowBox[{"ejecta", ":", " ", 
     RowBox[{
     "Artemeiva", " ", "2013", " ", "LPSC", " ", "abstract", " ", "TYCHO", 
      " ", "CRATER", " ", "EJECTA"}]}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"TempMel", "=", "1450.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"TempSod", "=", "1270.0"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"DensitySphere", "=", "3150.0"}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "Specific", " ", "heat", " ", "capacity", " ", "at", " ", "a", " ", 
     "fixed", " ", "pressure"}], ",", " ", 
    RowBox[{"unit", " ", "is", " ", "J", " ", 
     RowBox[{"kg", "^", 
      RowBox[{"-", "1"}]}], " ", 
     RowBox[{"K", "^", 
      RowBox[{"-", "1"}]}]}]}], "\n", "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"HeatCapacityConstPressure", "=", "1246.8"}], ";"}], 
  "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Stephan", " ", "Boltzmann", " ", "Constant"}], ",", " ", 
    RowBox[{"unit", " ", "is", " ", "J", " ", 
     RowBox[{"m", "^", 
      RowBox[{"-", "2"}]}], " ", 
     RowBox[{"s", "^", 
      RowBox[{"-", "1"}]}], " ", 
     RowBox[{"K", "^", 
      RowBox[{"-", "4"}]}]}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"StephanBoltzmanConst", "=", 
    RowBox[{"5.67", "*", 
     RowBox[{"10", "^", 
      RowBox[{"(", 
       RowBox[{"-", "8"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"Emissivity", " ", "of", " ", "sphere"}], ",", " ", "unitless"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Emissivity", "=", "0.7"}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Diameter", " ", "of", " ", "sphere"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"DiameterSphere", "=", "0.1"}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"Physical", " ", "constant", " ", "bundle"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"physicalconst1", "=", 
   RowBox[{"6.0", "*", "Emissivity", "*", 
    RowBox[{"StephanBoltzmanConst", "/", 
     RowBox[{"(", 
      RowBox[{
      "DensitySphere", "*", "HeatCapacityConstPressure", "*", 
       "DiameterSphere"}], ")"}]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"physicalconst2", "=", "TempEnv"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"TempIntegralPart", "=", 
   RowBox[{"Re", "[", 
    RowBox[{
     RowBox[{"2.0", "*", 
      RowBox[{"ArcTan", "[", 
       RowBox[{"TempSod", "/", "TempEnv"}], "]"}]}], "-", 
     RowBox[{"2.0", "*", 
      RowBox[{"ArcTan", "[", 
       RowBox[{"TempMel", "/", "TempEnv"}], "]"}]}], "-", 
     "\[IndentingNewLine]", 
     RowBox[{"Log", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"TempEnv", "-", "TempSod"}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{"TempEnv", "+", "TempSod"}], ")"}]}], "]"}], "+", 
     "\[IndentingNewLine]", 
     RowBox[{"Log", "[", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{"TempEnv", "-", "TempMel"}], ")"}], "/", 
       RowBox[{"(", 
        RowBox[{"TempEnv", "+", "TempMel"}], ")"}]}], "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"TimeCooling", "=", 
    RowBox[{
     RowBox[{"1.0", "/", 
      RowBox[{"(", 
       RowBox[{"physicalconst1", " ", "*", " ", "4.0", " ", "*", " ", 
        RowBox[{
         RowBox[{"(", "physicalconst2", ")"}], "^", "3"}]}], ")"}]}], "*", 
     "TempIntegralPart"}]}], ";"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"pos1", "=", 
     RowBox[{"Rg", " ", "*", " ", "0.54"}]}], ";"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"Print", "[", 
  RowBox[{
  "Df", ",", "\"\< \>\"", ",", "Rs", ",", "\"\< \>\"", ",", "Rg", ",", 
   "\"\< \>\"", ",", "pos1", ",", "\"\< \>\"", ",", 
   RowBox[{"pos1", "/", "Rg"}], ",", " ", "\"\< \>\"", ",", " ", 
   RowBox[{"vef", "[", "pos1", "]"}], ",", "\"\< \>\"", ",", 
   RowBox[{"Timeballistic", " ", "[", "pos1", "]"}], ",", "\"\< \>\"", ",", 
   "TimeCooling"}], "]"}]}], "Input",
 CellChangeTimes->{{3.6952381108921766`*^9, 3.6952382108738585`*^9}, {
   3.695238240935444*^9, 3.6952382576744585`*^9}, {3.695238315785203*^9, 
   3.6952384068279705`*^9}, {3.69523849414229*^9, 3.6952385560438833`*^9}, {
   3.695238608491756*^9, 3.6952386203791084`*^9}, {3.695238698598511*^9, 
   3.695238811387957*^9}, {3.6952388763939905`*^9, 3.695238963146703*^9}, {
   3.695246976255597*^9, 3.695247167481623*^9}, {3.6952472712378883`*^9, 
   3.695247296900053*^9}, {3.6952473308770704`*^9, 3.6952474265368834`*^9}, {
   3.6952474669879427`*^9, 3.695247474975194*^9}, {3.6952475339903727`*^9, 
   3.6952476818637204`*^9}, {3.6952478071169233`*^9, 3.695247857302445*^9}, {
   3.695247888658646*^9, 3.6952479563006797`*^9}, {3.6952481048136315`*^9, 
   3.6952482004266443`*^9}, {3.6952482842303815`*^9, 
   3.6952484440690064`*^9}, {3.6952484924761167`*^9, 
   3.6952484950969334`*^9}, {3.6952485470608664`*^9, 3.695248554611315*^9}, {
   3.6952486906129866`*^9, 3.695248714215938*^9}, {3.695248877221383*^9, 
   3.6952489629595327`*^9}, {3.695249006873814*^9, 3.695249008979828*^9}, {
   3.695249160222797*^9, 3.6952491717044706`*^9}, {3.695249236507286*^9, 
   3.69524923711569*^9}, {3.6952497863328104`*^9, 3.695249787690019*^9}, {
   3.695249820481429*^9, 3.6952498243970547`*^9}, {3.695249918418857*^9, 
   3.6952499330049505`*^9}, {3.695249997651765*^9, 3.6952501485359325`*^9}, {
   3.6952504763096333`*^9, 3.6952504764656343`*^9}, {3.6952511574255996`*^9, 
   3.69525118564618*^9}, {3.6952516667688646`*^9, 3.695251689311009*^9}, {
   3.695251828807103*^9, 3.6952518300863113`*^9}, {3.6952518773078136`*^9, 
   3.695251913344045*^9}, 3.6952519601443453`*^9, 3.695251992233751*^9, {
   3.695252024089155*^9, 3.6952521048196726`*^9}, {3.6952521659408646`*^9, 
   3.6952523176674366`*^9}, {3.6952524816868887`*^9, 3.695252481796089*^9}, {
   3.705754330605983*^9, 3.705754332103612*^9}, {3.7057548006742225`*^9, 
   3.7057548186145678`*^9}, {3.705754850782386*^9, 3.705754851453199*^9}, {
   3.705754959251272*^9, 3.705755024429325*^9}, {3.7057550721194425`*^9, 
   3.705755106408902*^9}, {3.7057653701068907`*^9, 3.705765392649324*^9}, {
   3.7057654279372025`*^9, 3.705765481227827*^9}, {3.7057655323500104`*^9, 
   3.7057655339568415`*^9}, {3.705766119326898*^9, 3.7057661912598815`*^9}, {
   3.7057662239425097`*^9, 3.7057662931766415`*^9}, {3.705766345999257*^9, 
   3.7057664565117826`*^9}, {3.706875478532056*^9, 3.7068755449728813`*^9}, {
   3.706875606016073*^9, 3.7068756064216757`*^9}, {3.7068756510379615`*^9, 
   3.7068756925030274`*^9}, {3.706879072972515*^9, 3.706879145232178*^9}, {
   3.706879802932394*^9, 3.706879803384797*^9}, {3.7068798861901283`*^9, 
   3.7068798864241295`*^9}, {3.7068799244727736`*^9, 3.706879928887602*^9}, {
   3.7068799678566513`*^9, 3.7068799689330587`*^9}, {3.70688083962144*^9, 
   3.7068808398866415`*^9}, {3.706880951864159*^9, 3.706880971613886*^9}, {
   3.706881033561883*^9, 3.7068810439983497`*^9}, {3.7068811542910566`*^9, 
   3.7068811652891273`*^9}, {3.7068812840370884`*^9, 3.706881285394297*^9}, {
   3.7068822117906356`*^9, 3.706882221103895*^9}, {3.7068823131288853`*^9, 
   3.706882316670108*^9}, {3.706882497054064*^9, 3.706882504354911*^9}, {
   3.7068825381447277`*^9, 3.706882540250741*^9}, {3.7068828601307917`*^9, 
   3.7068828611915984`*^9}, {3.7068828973526297`*^9, 
   3.7068829016270576`*^9}, {3.706883324623769*^9, 3.706883383888549*^9}, 
   3.7068834209231863`*^9, {3.7068835390939436`*^9, 3.706883539452746*^9}, {
   3.7068836296525245`*^9, 3.706883631009733*^9}, {3.7068936700145783`*^9, 
   3.706893694771937*^9}, {3.7068973706639366`*^9, 3.706897374111581*^9}, {
   3.7068975382561316`*^9, 3.7068975426865597`*^9}, {3.70706760480551*^9, 
   3.70706762979671*^9}, {3.70706769604991*^9, 3.70706774551751*^9}, {
   3.70706777585951*^9, 3.7070678444215097`*^9}, {3.70706790853751*^9, 
   3.70706791100231*^9}, {3.70706796408911*^9, 3.7070680762843103`*^9}, {
   3.70706812347431*^9, 3.70706826142511*^9}, {3.70706958859511*^9, 
   3.70706961439751*^9}, {3.70706965118231*^9, 3.70706968216391*^9}, {
   3.70707006019711*^9, 3.70707007116391*^9}, {3.70707017384311*^9, 
   3.70707017509111*^9}, {3.70707020992591*^9, 3.70707023550991*^9}, {
   3.70707036137071*^9, 3.70707042305311*^9}, {3.70707046665511*^9, 
   3.70707047790271*^9}, {3.7070705133147097`*^9, 3.70707059569831*^9}, {
   3.70707065383951*^9, 3.7070706784095097`*^9}, {3.7070707385631104`*^9, 
   3.70707093779071*^9}, {3.70707133270471*^9, 3.7070713586319103`*^9}, {
   3.70707139912951*^9, 3.7070713995507097`*^9}, {3.70707147404071*^9, 
   3.70707150586471*^9}, {3.710679729307892*^9, 3.7106797527078915`*^9}, {
   3.710680048842692*^9, 3.710680068451892*^9}, {3.710680175920292*^9, 
   3.7106801916918917`*^9}, {3.710680234747892*^9, 3.710680235683892*^9}, {
   3.710680304667092*^9, 3.710680304994692*^9}, {3.710680837984292*^9, 
   3.710680882428692*^9}, {3.7155315447696285`*^9, 3.715531551699028*^9}, {
   3.715531698553628*^9, 3.715531769344428*^9}, {3.715531820310828*^9, 
   3.7155318246516285`*^9}, {3.715531941336628*^9, 3.7155319427884283`*^9}, {
   3.715532060517628*^9, 3.715532090826828*^9}, {3.715532137571228*^9, 
   3.7155321398488283`*^9}, {3.715532830772828*^9, 3.715532850600428*^9}, {
   3.7156294653458605`*^9, 3.7156294870301385`*^9}, {3.715629523612608*^9, 
   3.715629524018213*^9}, {3.71579845430511*^9, 3.7157985038510275`*^9}, {
   3.7157985495593204`*^9, 3.7157985546137533`*^9}, {3.7157985924283953`*^9, 
   3.715798672300907*^9}, {3.715876159143181*^9, 3.7158761631835814`*^9}, {
   3.7158828383797812`*^9, 3.7158828674893813`*^9}, {3.715882904242981*^9, 
   3.715882908532981*^9}, 3.715883143718581*^9, {3.7160346080648336`*^9, 
   3.716034609968034*^9}, 3.7160346592952337`*^9, {3.7160347172492332`*^9, 
   3.7160347560308332`*^9}, {3.7160348074328337`*^9, 
   3.7160348633900337`*^9}, {3.7160349405476336`*^9, 
   3.7160349535268335`*^9}, {3.7160350093748336`*^9, 3.716035030746834*^9}, {
   3.716035192128834*^9, 3.7160352312380333`*^9}, {3.7160353242452335`*^9, 
   3.7160353518572335`*^9}, {3.7160354165192337`*^9, 
   3.7160354251616335`*^9}, {3.730898840790874*^9, 3.7308988488561773`*^9}, {
   3.730898891132719*^9, 3.7308989343764734`*^9}, {3.730902923690818*^9, 
   3.7309030759799705`*^9}, {3.730903159160237*^9, 3.730903243214114*^9}, {
   3.730903360184414*^9, 3.730903368608522*^9}, {3.73090342620446*^9, 
   3.730903426641266*^9}, {3.7309035778228045`*^9, 
   3.730903578103608*^9}},ExpressionUUID->"42fcca41-68e3-4358-9287-\
3a15a09535aa"],

Cell[BoxData[
 InterpretationBox[
  RowBox[{
  "98.17854670215006`", "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]",
    "38104.90978619579`", "\[InvisibleSpace]", "\<\" \"\>", 
   "\[InvisibleSpace]", "38569.23810788138`", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", 
   "15261.679622941754`", "\[InvisibleSpace]", "\<\" \"\>", 
   "\[InvisibleSpace]", "0.3956956468845343`", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", 
   "789.9655408000962`", "\[InvisibleSpace]", "\<\" \"\>", 
   "\[InvisibleSpace]", "688.7669430375416`", 
   "\[InvisibleSpace]", "\<\" \"\>", "\[InvisibleSpace]", 
   "88.0956896522931`"}],
  SequenceForm[
  98.17854670215006, " ", 38104.90978619579, " ", 38569.23810788138, " ", 
   15261.679622941754`, " ", 0.3956956468845343, " ", 789.9655408000962, " ", 
   688.7669430375416, " ", 88.0956896522931],
  Editable->False]], "Print",
 CellChangeTimes->{{3.7309030190392404`*^9, 3.730903078039197*^9}, {
   3.730903164245902*^9, 3.730903183324947*^9}, 3.730903220344221*^9, 
   3.730903370870551*^9, 3.730903429652105*^9, 
   3.730903581551252*^9},ExpressionUUID->"3df6e255-d9ee-4db9-825a-\
81d9162d3c4f"]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.70706779535951*^9, 
  3.70706780275391*^9}},ExpressionUUID->"0aa5f1fa-9f4a-4dd5-90d3-\
b4af8e408de2"],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.705765437281782*^9, 
  3.7057654404174423`*^9}},ExpressionUUID->"8f8e12ec-6bd8-4ca2-ae91-\
a63f75e463f5"]
},
WindowSize->{789, 736},
WindowMargins->{{6, Automatic}, {Automatic, 0}},
FrontEndVersion->"11.1 for Microsoft Windows (64-bit) (April 18, 2017)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 20574, 522, 1456, "Input", "ExpressionUUID" -> \
"42fcca41-68e3-4358-9287-3a15a09535aa"],
Cell[21157, 546, 1167, 23, 22, "Print", "ExpressionUUID" -> \
"3df6e255-d9ee-4db9-825a-81d9162d3c4f"]
}, Open  ]],
Cell[22339, 572, 150, 3, 30, "Input", "ExpressionUUID" -> \
"0aa5f1fa-9f4a-4dd5-90d3-b4af8e408de2"],
Cell[22492, 577, 154, 3, 30, "Input", "ExpressionUUID" -> \
"8f8e12ec-6bd8-4ca2-ae91-a63f75e463f5"]
}
]
*)

